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Abstract 

It is shown that if a complete set of mutually commuting operators is formed by constants of 
motion, then, up to a factor that only depends on the time, each common eigenfunction of such 
operators is a solution of the Schrodinger equation. In particular, the operators representing the 
initial values of the Cartesian coordinates of a particle are constants of motion that commute 
with each other and from their common eigenfunction one readily obtains the Green function. 

PACS numbers: 03.65.-w 


1 Introduction 

In the standard procedure to solve the Schrodinger equation for a time-independent Hamiltonian 
one makes use of the method of separation of variables, or one starts looking for a set of mutually 
commuting operators that also commute with the Hamiltonian. In fact, the separable solutions 
of the Schrodinger equation are common eigenfunctions of a set of mutually commuting operators, 
with the separation constants being the eigenvalues of such operators. Usually, these operators do 
not depend explicitly on the time and, since they commute with the Hamiltonian, are constants of 
motion. 

As we shall show below, there is no reason to restrict ourselves to time-independent operators; we 
can find solutions of the Schrodinger equation that are common eigenfunctions of a set of, possibly 
time-dependent, mutually commuting operators that are constants of motion (and, therefore, may 
not commute with the Hamiltonian). This method is analogous to that given by the Liouville 
theorem on the solutions of the Hamilton-Jacobi equation (see, e.g.. Ref. [I]). 

In Section 2 we establish the basic results of this paper and in Section 3 we present several 
examples, exposing the advantages of the method. We show that making use of the operators that 
represent the initial position of a particle one readily obtains the Green function of the corresponding 
Hamiltonian. 


2 Basic results 

Let A be a, possibly time-dependent, Hermitean operator that is a constant of motion, i.e., 

m^ + lA,N] = o, ( 1 ) 

where H is the Hamiltonian of the system under consideration. If b is an eigenfunction of A with 
eigenvalue A, 

Ab = Ab (2) 

and il> satisfies the Schrodinger equation 

= Hi;, (3) 


1 


then A is a constant. Indeed, taking the derivative with respect to the time of both sides of m , we 
obtain 


_dA 


dt 

and, making use of Eqs. m and dni), 




ih-— \ ^1) A[ ih-— = + Xih-;— 


dt 


AX 




dt 


dt 


dA 


that is, 


— [A, H]4’ + AHijj = + Ai7V') 


HAip = ifi,—i/) + XHtp 
at 


and using Eq. ([2|) again, it follows that dA/dt = 0. 

Conversely, assuming that A is a constant of motion and that its eigenvalues are constant, from 
© we have 




^di^X _ di, 


i?i-^ J Ip -\- A\ ) = Xih 


dt 


dt 


or, by virtue of o. 


which amounts to 


or 


-[A H]'tp + A{ = Ah^, 

-AHiP + HXip + A ( ih^) = XiH^, 


dt 


dt 


.^dip 


A ( '\h-^ — Hip ) = A ( ~ Hip ) . 


Alp 


dt 

If the spectrum of A is non-degenerate, it follows that 

dip 

-Hip = lip), 

where /i is some complex-valued function of t only. Then, if F{t) is a solution of 

t^F -k ih— = 0, 

we find that Fip satisfies the Schrddinger equation. In fact. 


(4) 


d{Fp,) 


.^dF\ 


= F{Hip + ^i'tp)+\ ih—j Ip 

dt ) 


= iJ(F^)-I-( -I-ili^ ] 


(Note that, if ip is an eigenfunction of A with eigenvalue A, then, if F{t) is a function of t only, 
F{t)ip is also an eigenfunction of A with eigenvalue A.) When the spectrum is degenerate, the proof 
is similar, considering linear combinations of the eigenfunctions of A with the same eigenvalue A and, 
in place of ([4]), one would have to solve a linear system of first-order differential equations, which 
may be as complicated as the original Schrddinger equation. However, if we find a complete set of 
mutually commuting operators that are constants of motion (complete in the sense that, up to a 
factor, there is only one common eigenvector of these operators for a given set of the eigenvalues), 
the common eigenfunctions of such a set of operators can be chosen in such a way that they are 
solutions of the Schrddinger equation (assuming that the eigenvalues are constant). 

This result is analogous to the Liouville theorem of the Hamiltonian mechanics, according to 
which if we have a set of constants of motion such that their Poisson brackets are equal to zero, they 
can be used to find complete solutions of the Hamilton-Jacobi equation. 


3 Examples 

In this section we give several examples, finding solutions of the Schrddinger equation starting from 
the eigenfunctions of constants of motion that depend explicitly on the time. 
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3.1 The one-dimensional harmonic oscillator 


One can readily verify that the operator 


Xq = X cos uit -sin ujt 

muj 


is a constant of motion if the Hamiltonian is given by 


H = 


P 


2m 2 

The eigenvalue equation XQ^p = xgip amounts to the first-order differential equation 

h sin cot d 


X cos uit — 


imuj dx 


= Xotp 


whose solution is readily found to be 

'll) = F{t) exp 


2hsinujt 


{x^ cos wt — 2xi)x') 


(5) 

( 6 ) 


(7) 


where Fit) is a function of t only. 

Assuming that xq is a constant, one finds that the wavefunction ([7]) satisfies the Schrodinger 
equation, ihdip/dt = Hip, if and only if [cf. Eq. (|4])] 


dlnF 

dt 


uj cos uit muP'Xif' 
2sina;t 2ihsiv? ujt' 


hence 


F = 


C 


\/sin uit 


exp 


imujXff cosuit 
2hsinuit 


where C is a constant, and substituting this expression into © we obtain 


V'xo (x, t) = ^ exp j Ux^ + cos uit - 2xqx\ \ . (8) 

Vsmwt 12h sin Lot ^ ■'J 

As usual, the wavefunctions ipxo are orthogonal for different values of the eigenvalue xq. Note also 
that the functions ([5|) are not separable. 

In the present example, F[ is time-independent and, therefore, there is a set of stationary states, 
(^„(x) exp(—iE„t/?i), where the are the eigenvalues of FI. Hence, the wavefunction (|S]) must be 
expressible in the form 

lpxo{x,t) = y^ Cn{xo)4ln{x) exp 
n 

and the symmetry of the right-hand side of (|H1) under the interchange of x and xq implies that 



Ipxoixf ) = '^(t)n{xo)(t)n{x) exp ( - 


iEr,.t 


i.e., with the appropriate value of the constant C, ip^oixf) is the time-dependent Green function of 
the one-dimensional harmonic oscillator (cf., for instance, Refs. [2l|3l|4]). This was to be expected 
since the operator Xq is, in the Heisenberg picture, the position operator at t = 0 and, therefore, 
■0X0 corresponds to the state of the particle with a well defined value (equal to Xq) of the position 
at t = 0. The value of the normalization constant C can be determined by the condition 

/ CO 

'ij}o{x,t) dx = I, (9) 

-oo 

which gives 

C=,[^. 

Y 27ri?i 

As is well known [3] , from the Green function (jS]) one can obtain the expression of the stationary 
states (j)n{x) exp{—iEnt/h) in terms of the Hermite polynomials. 
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3.2 A time-dependent Hamiltonian 

Now we shall consider the time-dependent Hamiltonian 

o2 


H = ^ -fete, 

2m 


where fc is a constant. The operator 


Xn = X- 


tp kt^ 


( 10 ) 


( 11 ) 


m 3m 

is a constant of motion and, in the Heisenberg picture, corresponds to the position of the particle at 
t = 0. The eigenfunctions of this operator are determined by the first-order differential equation 

ht d kt^\ ^ , 

^ ]V = XQip, 

im ax 3m ) 

where xq is the eigenvalue of Xq. The solutions of this equation have the form 

kt^x 


Ip = F{t) exp 


\m X 


ht 


2 3m 


where F{t) is a function of t only. Substituting this expression into the Schrodinger equation, 
assuming that xq is constant, we obtain the equation 

dlnF 


dt 


thus. 


C 


1 m 
2t 2iht‘^ 


„ dkt^xo 
xo + 


3m 


9m'^ J 




(x - Xo) + 


2 2kF{x — Xo) 4k^t 


2x61 


3m 


45m^ 


( 12 ) 


where C is a constant. With the appropriate value of C (C = •\/m/27ri?i), (IT^ is the Green function 
for the Hamiltonian (flOt . 

It may be remarked that, since the Hamiltonian (1101) depends explicitly on the time, the corre¬ 
sponding Schrodinger equation cannot be solved by separation of variables. 


3.3 Another standard example 

Another example usually considered in connection with the Green functions is that of a particle in 
a uniform field. If 

H = eEx, (13) 

2m 

where e and E are constants, one readily finds that 

tp eEt^ 

Xo = x-^ + ^— (14) 

m 2m 

is a constant of motion that corresponds to the position of the particle at f = 0. Thus, the eigenvalue 
equation Xoip = xoip is equivalent to 


ht d eEt^ 
im dx 2m 


with the solution 


Ip = E(t) exp 


im 

2ht 


X — 2xox + 


tp = Xoip, 
eEt^x 


where E[t) is a function of t only. Substituting this last expression into the Schrodinger equation 
one obtains the condition 


and, therefore. 


dlnF 

dt 


C 


m 


2t 2ife2 


2 eEt^xo e^E'^t^ 
Xo - 


V'rro (2;. 0 ^ exp <^ — 


Vt 


2ht 


[x - xo)^ -b 


m 4m^ 

eEt^{x + xo) e^E'^t'^'^ 
m 12m? 


where C is a constant. 
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3.4 A two-dimensional system 


The Hamiltonian 




+ 



(15) 


corresponds to a charged particle of mass m and electric charge e in a uniform magnetic field B. 
The operators 


^0 


>1) 


I, ,1 1 1 , 

-(1 + cosLutjx — -ysmujt -pjjSinwt H-(1 — cosujtjpy, 

2 2 moj m(jj 

1. N 1 1 1 / 

- (1 + cos uit)y + —X sin cot - py sin tot -(1 — cos ujt)px , 

2 2 moj moj 


(16) 

(17) 


where w = eB jmc^ are two constants of motion that commute with each other. The operators p6l) 
and (O correspond to the Cartesian coordinates of the particle at t = 0. (It may be remarked that, 
for any system, the operators that represent the initial values of the Cartesian coordinates, or of the 
Cartesian components of the momentum, are constants of motion that commute with each other.) 

Since the operators m and (II3 are linear in the momenta, one readily finds their eigenfunctions. 
In fact, the eigenfunctions of Xq and Yq with eigenvalues xq and y^, respectively, are 


V' 


F{t) exp 
F{t) exp 


Ah{l — cos wt) 


4S 


[(^ 


[{x'^ + y^ — 2xox — 2yoy) sinuit + 2{xoy — yox)il — cos cat)] j 
- 2xox - 2yoy) cot + 2{xoy - yox)] i , (18) 


where F{t) is a function of t only. Substituting (1181) into the Schrodinger equation we obtain the 
single condition 


dlnT" 

dt 


UJ , imW^ , 2 2^ 2 1 

--cot^wt-§^(^0 +2^0 2^^’ 


therefore. 




C CSC ^ojt exp 


imw 


{[(a: - xo)^ + (j/ - yo)^] cot ^ujt + 2{xoy - yox)} , 


where C is a constant. 


(19) 


4 Concluding remarks 

The procedure followed here to find the Green functions highly contrasts with the methods usually 
employed in the literature (see, e.g.. Refs. mMm)- For instance, the fact that the operators Xq, 
defined by Eqs. ®, dni), and m, are linear in p implies that, in order to find the Green func¬ 
tions, one only has to solve two first-order ordinary differential equations. Even though the results 
presented in Section 2 also apply in the usual case of time-independent operators that commute 
with the Hamiltonian, the difference is that there are not many useful time-independent constants 
of motion linear in the momentum. 
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